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Abstract 



Working within the framework of the minimal supergravity model with 
gauge coupling unification and radiative electroweak symmetry breaking, we 
evaluate the cosmological relic density from lightest neutralinos produced in 
the early universe. Our numerical calculation is distinct in that it involves 
direct evaluation of neutralino annihilation cross sections using helicity am- 
plitude techniques, and thus avoids the usual expansion as a power series 
in terms of neutralino velocity. Thus, our calculation includes relativistic 
Boltzmann averaging, neutralino annihilation threshold effects, and proper 
treatment of integration over Breit-Wigner poles. We map out regions of 
parameter space that give rise to interesting cosmological dark matter relic 
densities. We compare these regions with recent calculations of the reach for 
supersymmetry by LEP2 and Tevatron Main Injector era experiments. The 
cosmologically favored regions overlap considerably with the regions where 
large trilepton signals are expected to occur at the Tevatron. The CERN 
LHC pp collider can make a thorough exploration of the cosmologically fa- 
vored region via gluino and squark searches. In addition, over most of the 
favored region, sleptons ought to be light enough to be detectable at both 
LHC and at a y/s = 500 GeV e+e" collider. 
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I. INTRODUCTION 



The Minimal Supersymmetric Standard Model [TJ (MSSM) is one of the leading can- 
didate theories for physics beyond the Standard Model (SM). The MSSM is a globally 
supersymmetric version of the SM, where supersymmetry breaking is implemented by the 
explicit introduction of soft-supersymmetry breaking terms. The MSSM is minimal in the 
sense that the fewest number of additional new particles and interactions are incorporated 
which are consistent with phenomenology. In particular, possible baryon (B) and lepton 
(L) number violating interactions are excluded from the superpotential (the presence of both 
B and L violating interactions can lead to catastrophic proton decay rates). As a result, 
there exists a conserved .R-parity, where the multiplicative quantum number R = +1 for 
ordinary particles, and R — — 1 for superpartners. A consequence of .R-parity conservation 
is that the lightest supersymmetric particle (LSP) is absolutely stable. Theoretical prejudice 
coupled with experimental constraints strongly favor a color and charge neutral LSP ||. In 
addition, in the MSSM, the LSP is strongly favored to be the massive, weakly interacting 
lightest neutralino Z\ fj|||. Z\&, if they exist, would have been abundantly produced in the 
early universe; if so, then relic neutralinos could well make up the bulk of the dark matter 
in the universe today |].[|. 

In this paper, we restrict ourselves to the framework of the low energy effective La- 
grangian which is expected to result from, for instance, supergravity grand unified models 
0. In these models, it is assumed that supersymmetry is broken in a hidden sector of 
the theory. Supersymmetry breaking is then communicated to the observable sector via 
gravitational interactions, leading to a common mass mo for all scalar particles, a com- 
mon mass mi/2 for all gauginos, a common trilinear coupling A , and a bilinear coupling 
B . These soft supersymmetry breaking parameters are induced at energy scales at or be- 
yond the unification scales, but with (theoretically motivated) values typically in the range 
100 — 1000 GeV. The resulting theory, the MSSM with universal soft breaking terms, is 
then regarded as an effective theory with Lagrangian parameters renormalized at an ultra- 
high scale Mx ~ Mqut — Mpi anc k, and valid only below this scale. The corresponding 
weak scale sparticle couplings and masses can then be calculated by evolving 26 parame- 
ters via renormalization group equations || from the unification scale to the weak scale. 
An elegant by-product |5| of this mechanism is that one of the Higgs boson mass squared 
terms is driven negative, resulting in a breakdown of electroweak symmetry. The radiative 
electroweak symmetry breaking constraint allows one to essentially eliminate B in favor of 
tan/3 (the ratio of Higgs field vev's), and to calculate the magnitude of the superpotential 
Higgs mixing term p in terms of Mz (where we actually minimize the full one-loop effective 
potential). The model is then specified by only four SUSY parameters (in addition to SM 
masses and couplings). A hybrid set consisting of the common GUT scale scalar mass mo, 
common gaugino mass m^, common SUSY-breaking trilinear coupling Aq, along with the 
weak scale value of tan (3 proves to be a convenient choice. In addition, the sign of p must 
be stipulated. These parameters fix the weak scale masses and couplings of all the sparticles 

The matter density of the universe p is usually parametrized fHH in terms of Q = p/p c , 
where p c = 3Hq/8tiGn — 1-88 x 10~ 29 /i 2 g/cm 3 , and h, the Hubble scaling constant, is 
related to the Hubble constant H by H = lOO/i km/sec/Mpc. Here h parametrizes our 
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ignorance of the true value of Hq, so that 0.5 ^ h ^ 0.8. Measurements of galactic rotation 
curves suggest Q ~ 0.03 — 0.1, compared to a luminous matter density of fli um . ~ 0.01. 
Galactic clustering and galactic flows suggest even larger values of Q ~ 0.2 — 1. Finally, the 
theoretically attractive inflationary cosmological models require a flat universe with Q = 1. 
Meanwhile, estimates of the baryonic contribution to the matter density of the universe 
from Big-Bang nucleosynthesis suggest that flbaryonic ~ 0.01 — 0.1. These analyses and 
estimates suggest that the bulk of matter in the universe is (non-baryonic) dark matter. 
Finally, analyses of structure formation in the universe in light of the COBE measurements 
of anisotropies in the cosmic microwave background radiation suggest that dark matter may 
be made of ~ 60% cold-dark matter (weakly interacting massive particles or WIMPS, such 
as the lightest neutralino Zi), ~ 30% hot dark matter (such as relic neutrinos), and ~ 10% 
baryonic matter. This is the so-called "mixed dark matter " scenario. 

The central idea || behind relic density calculations is that in the very early universe, 
neutralinos were being created and annihilated, but that they were in a state of thermal 
equilibrium with the cosmic soup. As the universe expanded and cooled, temperatures 
dropped low enough that neutralinos could no longer be produced (T < although 
they could still annihilate with one another, at a rate governed by the total neutralino pair 
annihilation cross section, and the neutralino number density. Ultimately, as the universe 
expanded further, the expansion rate outstripped the annihilation rate, thus freezing out the 
remaining neutralino population of the universe, and locking in a neutralino relic density. 
Our goal in this paper is to carry out estimates of the neutralino relic density expected 
from the minimal SUGRA model. One solid constraint on supersymmetric models with 
relic dark matter particles comes from the age of the universe, which ought to be greater 
than 10 (15) Gyrs; this implies Qh 2 < 1 (0.25). Thus, models with too large a relic density 
would yield too young of a universe, in violation at least with the age of the oldest stars in 
globular clusters. Furthermore, models with Qh 2 < 0.025 would not even be able to account 
for the dark matter needed to explain galactic rotation: such models would be considered 
cosmologically uninteresting. Models with intermediate values of 0.025 ^ Qh 2 < 1 are 
considered cosmologically interesting, as they might explain galactic rotation and clustering, 
or might even make up the matter density needed for inflationary cosmology, given a cold- 
dark matter (CDM: Qh 2 ~ 0.25 — 0.64) or mixed hot/cold dark matter scenario (MDM: 
ft/i 2 ~ .15- .4). 



Following the procedures outlined by Lee and Weinberg ||11|| , many groups have calcu- 
lated the relic neutralino abundance |T2Trn|. Early works involved calculating the most 
important neutralino annihilation channels, usually assuming the LSP was a photino. Later 
studies included various improvements, including more annihilation channels, more general 
neutralino mixings, and more realistic supersymmetric particle spectra. A common thread 
amongst many papers was the calculation of the Boltzmann-averaged quantity a x v using 
a power series expansion in velocity. Such an approach was shown to be inaccurate when 
relativistic effects were important, when annihilation proceeded through s-channel poles, 
when threshold effects were important, or when co-annihilation processes occured pU| , |2"T|j . 
Many recent calculations have included some or all of these effects. 

We have several goals in mind for the present paper. 

• We wish to present reliable calculations for the neutralino relic density in supersym- 
metric models. To this end we evaluate all 2 — > 2 neutralino annihilation diagram 
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amplitudes numerically as complex numbers, without approximation. We perform 



Boltzmann averaging using the Gondolo-Gelmini formalism ||21|| . This takes into ac- 
count relativistic thermal averaging, while our numerical helicity amplitude technique 
avoids the usual uncertainties inherent in the velocity expansion, so that Breit-Wigner 
poles and threshold effects are fully accounted for. Co- annihilation can occur when the 
two lightest superpartners are very close in mass- this situation rarely occurs within 
the SUGRA framework adopted in this paper, and hence we ignore it. 

We present results in the well-motivated SUGRA framework, which includes gauge 
coupling unification, Higgs mass radiative corrections ||22|| , and radiative electroweak 



symmetry breaking using the one- loop effective potential [ f23|j . 

Our results for the relic density calculation can be directly related to recent calculations 
for various supersymmetry signals expected at the LEP2 [p4,25], Tevatron |2~6H2~9fl and 



LHC colliders []30| , [31}| . In particular, relic density calculations have a preference for 
light sleptons. Such light sleptons may well be observable at LEP2 or LHC colliders, 
and yield enhanced rates for W\Z 2 — > 3£ states at the Tevatron collider |32| . 



To accomodate these goals, we present in Sec. II various details of our relic density cal- 
culation, including those peculiar to the present approach. In Sec. Ill, we present numerical 
results for our relic density calculations in the tuq vs. plane of the minimal SUGRA 
model. In Sec. IV, these results are explicitly compared to expectations for minimal SUGRA 
at various collider experiments, as worked out in a series of previous papers. Finally, in Sec. 
V we present an overview and some conclusions. 



II. CALCULATIONAL DETAILS 

We begin our determination of the neutralino relic density from minimal supergravity 
by selecting a point in the SUGRA parameter space 

m o, m i/2, tan/?, A and sign(fi), (2.1) 

where in addition we take the top quark mass m t = 170 GeV. The 26 renormalization group 
equations are iteratively run between the weak scale and the GUT scale, which is defined 
as the point where the U(l), SU(2) and SU(3) gauge couplings unify, and is typically 
M x ~ 2 x 10 16 GeV. We use 2-loop RGE equations for gauge and Yukawa couplings (with 
SUSY particle threshold effects), but only 1-loop equations for the running of the various 
soft breaking terms. The 1-loop effective Higgs potential is minimized to enforce radiative 
electroweak symmetry breaking. Our procedure has been described in more detail in Ref. 
33fl , and has been incorporated into the event generator ISA JET ||34|| . At this point, a 



correlated sparticle mass spectrum and couplings emerge from our input point in SUGRA 
parameter space. 

The next step in our computation, after obtaining the superparticle spectrum, is to 
evaluate the neutralino relic density by solving the Boltzmann equation as formulated for a 
Friedmann- Robertson- Walker cosmology |5[]. Central to the evaluation of the relic density 
is the computation of the fully relativistic, thermally averaged neutralino annihilation cross 
section times velocity, defined as 
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< av Mol > (T) - j e - El/Te -E 2/ T d s pid s p2 > ( 2 - 2 ) 

where pi and P2 (-^2) are the momentum and energy of the two colliding particles in 
the cosmic, co-moving frame of reference, and T is the temperature. The above expression 
has been reduced to a one-dimensional integral by Gondolo and Gelmini [21|, which yields 



< crvuol > O) = 4x ^2/-i\ J da°(a)a 2 (a 2 ~ 4 )^i(~)' ( 2 - 3 ) 



where x = a = \fs is the subprocess energy, and Ki are modified Bessel functions 

m z 1 m z 1 

of order i. 

We evaluate the neutralino annihilation cross section for Z\Z\ — > fif 2 as 

<M« = on ,1 2 f ^l^cosfl, (2.4) 
6Z7TS \ 2<s,m~ ,m~ ) 

where | 2 is the spin summed and averaged squared matrix element. Our calculation of 
the relic density is distinct in that we evaluate M. for all Feynman diagrams listed in Table 
1 as complex numbers, using the HELAS JHJ helicity amplitude subroutine package. Thus, 
our approach avoids the usual uncertainties associated with the expansion of cross section 
in terms of a power series in velocity. The integration over cos# is performed numerically 
using Gaussian quadratures. 

To evaluate the neutralino relic density, the freeze out temperature xp is needed. The 
standard procedure here to iteratively solve the freeze out relation 



\ m z / 45 in 

F = lo § TT~3 V o r< < aVM ° l >x f X F > ( 2 - 5 ) 
L An- 3 y Zg*GN J 

by starting with a trial value xp = In the above, g* is the effective number of degrees of 
freedom at T = Tp (y/g^ ~ 9), and Gn is Newton's constant. 
Finally, the relic density can be calculated from 

^ 2 = 1, (2-6) 

8.0992 x 10" 47 GeV 4 V ; 



where 



T 



1 



p(T ) 1.66 x {-^fT^-— — -. (2.7) 

Mpi T 7 1 Jo < ovuoi > dx 

To evaluate the integral in the above expression, we expand the modified Bessel functions 
in Eq. 2.3 as power series in x, and then integrate over x. The result is 

r x F I 



rxp I roo 3 

/ < avMoi > dx = — — / daa(a)a^ (a 2 — 4)F(a), (2.8) 
Jo on J2 



where 
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2 285 45 15 
3^^2~ ~ 32^ ~ 18^ X 

{ g-2 
e ^ 

In the above, virtually all the contribution to the integral comes from x < 2.5. We integrate 
the above expression numerically with Gaussian quadratures, taking care to scan finely the 
regions with a Breit-Wigner pole. In the region of a pole, the domain of integration must be 
broken into very tiny intervals, and obtaining convergence for a single point in parameter 
space can take up to several hours of CPU time on a DEC ALPHA. 



x 2 F -2{ 



+ 2v^(a- 2)5(1 -Erf( 



xf 
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III. RESULTS FROM RELIC DENSITY CALCULATION 

Our first numerical results for the relic density from minimal SUGRA models are given 
in Fig. 1, where we plot contours of the neutralino relic density flh 2 in the mo vs. m\/2 
parameter plane, where we take A = 0, tan/3 = 2, /i < and m t = 170 GeV. Changes in 
the Aq parameter mainly affect 3rd generation sparticle masses, and consequently result in 
only small changes in the relic density. The regions labelled TH are excluded by theoretical 
considerations: either there is a charged or colored LSP (or the v is LSP), or the radiative 
electroweak symmetry breaking constraint breaks down. The region labelled by EX corre- 
sponds to parameter space already excluded by SUSY searches at LEP or Fermilab Tevatron 
experiments H55f| . 



In almost all of the plane, we find Qh 2 > 0.025, i.e. large enough to explain the galactic 
rotation curves. However, the region to the right of the Qh 2 = 1 contour is certainly excluded 
in that the age of the universe would be younger than 10 Gyrs. Meanwhile, a dominantly 
CDM inflationary universe would lie in between the flh 2 = 0.25 — 0.75 contours. The 
COBE favored MDM inflationary universe would lie between the Qh 2 = 0.15 — 0.4 contours. 
For this latter favored region, mi/ 2 is bounded by m\/2 ^ 400 GeV (corresponding to 
rrig < 1000 GeV), and m < 150 GeV, unless the gluino is very light {m g ~ 300 GeV). (For 
comparison, various SUSY particle mass contours for the same parameter choices are listed 
in Refs. p3| , p8| , p0| .) We find in general that large values of m ^ 350 GeV (corresponding to 
m~ t > 250 GeV) yield too young a universe (due to suppression of t-channel slepton exchange 
diagrams), except for the two narrow corridors in the lower right region of the figure. The 
upper of the two corridors corresponds to neutralino annihilation through the Z pole, so the 
relic density is largely reduced by Z mediated s-channel annihilation diagrams. The lower 
of the two corridors corresponds to annihilation through an s-channel light Higgs h pole- in 
this case, the relic density falls rapidly to values even below Qh 2 ~ 0.025. 

A qualitative feel for the relative importance of different annihilation channels can be 
gleaned from Fig. 2. Here we plot for itlq fixed at 200 GeV, as a function of rrii/2, the 
thermally averaged annihilation cross section times velocity, integrated over temperature, 
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which enters into the relic density calculation (Eq. 2.8). Larger cross sections correspond 
to smaller relic densities. As m^j increases, the first pole we come to is annihilation via 
s-channel h, where Z 1 Z 1 — > bb dominates. In these plots, scales with m^, and at 
the Higgs pole in this plot (on the edge of exclusion by LEP Higgs search experiments), 
m-^ — 30 GeV and m^ — ^ GeV. As one moves to higher m.x/2, annihilation through 
the Z pole is reached, which is dominated by Z X Z\ — ► dd, ss, and bb. For values of mx/2 
away from poles, annihilation via t-channel slepton and sneutrino exchange dominates. For 
even higher mx/2 values, annihilation into channels such as hh, Zh, WW and ZZ open up, 
but never dominate for the parameter choices in this plot. Annihilation into other channels 
such as HA, AA, HH amd H + H~ are included in our calculation, but unimportant given 
our SUGRA sparticle mass spectrum, which yields very large masses for Higgs bosons other 
than h. The onset of the Z\Zx — > tt can be detected in the Z\Z\ — > T,UiUi curve around 
mi/2 ~ 400 GeV. 

If we plot the relic density for the same parameter choices, but flip the sign of n, so that 
H > 0, then we obtain the results of Fig. 3. The relic density contours in this case are 
similar to those of Fig. 1 for large values of miw, where annihilation dominantly occurs via 
slepton exchange. The kink in the contours is due to the onset of the Z\Z\ — ► tt channel. In 
this case, annihilation through t-channel t\ exchange makes a large contribution to the total 
annihilation cross section. For smaller values of m^, in contrast to Fig. 1, we find only one 
corridor extending to large m where the relic density drops to cosmologically un-interesting 
values. In this case, the Z and h poles very nearly overlap for ~ 46 GeV. This can 
be seen in more detail in Fig. 4, where again we show the thermally averaged cross section 
versus mx/2, for m Q = 200 GeV. 

Finally, we show again the neutralino relic density flh 2 in the mo vs. mx/2 plane for the 
same parameter choices as Fig. 1, except now we take a large value of tan/? = 10. For this 
case, we note the rather broad band at mx/2 ~ 100 — 140 GeV, where Qh 2 < 0.025- too 
low to explain even the galactic rotation curves, and due again to annihilation through the 
s-channel graphs. In fact, inflationary models, which require Qh 2 > 0.15, are only allowed if 
mx/2 > 150 GeV, corresponding to rrig > 400 GeV. In this plot, there is a significant region 
extending to large values of mo, corresponding to large rriq and large m^, for m\j2 ~ 150 — 190 
GeV. The contributing thermally averaged subprocess cross sections are again shown in 
Fig. 6 for m = 200 GeV. In this plot, the Z pole annihilation channel occurs at m^j — 110 
GeV, followed by the Higgs pole at mi/2 — 130 GeV. The rough overlap of these two pole 
contributions leads to the single broad corridor of low flh 2 shown in Fig. 5. 



IV. IMPLICATIONS FOR SUSY SEARCHES AT COLLIDERS 



Recently, various papers have been written on the prospects for supersymmetry at the 
LEP2 e+e" collider |2§|25|, the Fermilab Tevatron pp collider and the CERN 



LHC pp collider p0|,31| . Our objective in this section is to assess the prospects for discovery 



of SUGRA at hadron and e + e~ colliders, given the additional constraints from requiring a 
reasonable value for the neutralino relic density. We mainly focus on the collider results of 
Refs. [^,^,^,^,^T|,^] , since they were performed in a consistent framework, in the same 
m vs. mx/2 plane. 
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In Fig. 7, we again show the neutralino relic density contours in the mo vs. pa- 
rameter plane, for the same parameter choices as in Fig. 1. In addition, we have added 
on contours for SUSY discovery at various colliders. Supersymmetric particles ought to be 
discoverable at LEP2 operating at y/s = 190 GeV, with integrated luminosity / Cdt = 300 
fb _1 below the contour labelled LEP2 [Z5[]. The lower-left bulge in the LEP2 contour is 
where sleptons ought to be detectable, while beneath the contour running along 777,1/2 — 100 
GeV (which runs through the neutralino Z-pole annihilation region), charginos ought to 
be detectable. By comparing, we see that the region accessible by LEP2 generally has 
Qh 2 < 0.15 i.e. not the most cosmologically favored region, but with enough dark matter to 
explain galactic rotation. However, the contour labelled with LEP2-Higgs shows the reach 
of LEP2 for the light SUSY Higgs boson h, which is just below 777,1/2 ~ 400 GeV. This region 
completely encloses the favored MDM region. The implication is that if MDM explains dark 
matter in the universe, and if tan/5 is small and /i < 0, then LEP2 ought to discover at least 
the light SUSY Higgs boson. 

The dashed contour labelled Tevatron is a composite of the reach of Tevatron Main 
Injector era (y/s = 2 TeV; J Cdt = 1 fb" 1 integrated luminosity) experiments for multi- 
jet+^ T events [p6] , and mainly, for W\Z 2 — > 3£ events p8fl . We see that the largest reach 



by Tevatron experiments occurs exactly in the cosmologically favored MDM region, and can 
reach to myj ~ 160 GeV, corresponding to m~ g ~ 440 GeV. This is no accident: a reasonable 
neutralino annihilation cross section generally requires m~ t < 200 GeV; these lighter sleptons 
give rise to enhanced leptonic decay of neutralinos, leading to large rates for W\Z 2 — > 3i 
events. Since lower values of 777,1/2 are preferred by fine-tuning arguments [ fjfi[| , there is a 
good chance Tevatron experiments could discover SUSY via 3i events if nature chose this 
parameter set. 

We also compare the results of Fig. 7 with expectations for supersymmetry at the CERN 
LHC collider. Of course, LHC experiments can cover the whole parameter plane up to 
mi/2 ~ 600 — 800 GeV with only / Cdt = 10 fb _1 of integrated luminosity, at y/s = 14 TeV, 
via searches for multi-jet+^ T events from gluino and squark cascade decays fl30l , so discovery 
of SUSY would be no problem. We also plot in Fig. 7 the contour beneath which sleptons 
ought to be visible at LHC |30"1 , |31~1| . We see that the cosmologically favored MDM region falls 
almost entirely within the slepton discovery region, so that if the MDM scenario is correct, 
then LHC has a very high probability to discover a slepton. Since sleptons are relatively 
light, LHC experiments ought as well to be sensitive to W\Z 2 — > 3C, events over much, but 
not all, of the favored MDM region [^,^1|]. (In some of the favored region, Z 2 — > vv or Z\h, 
thus spoiling the signal.) Finally, since < 250 GeV in the MDM scenario, sleptons would 
then likely be visible at a linear e + e~ collider operating at y/s = 500 GeV. 

In Fig. 8, we show the same relic density contours as in Fig. 3 (tan/3 = 2 , /x > 0), and 
compare again with expectations for colliders. In this case, we see the LEP2 contour again 
lies in a region of Qh 2 < 0.15, although it does encompass the cosmologically interesting 
region around (m Q ,m 1 / 2 ) ~ (100, 110). The LEP2 Higgs contour in this case lies at ~ 
170 GeV, and thus covers only a portion of the MDM favored region. Thus, if the MDM 
scenario is correct, and tan /3 is small, minimal SUGRA sparticles or light Higgs boson might 
still not be accessible at LEP2. We also plot the contour due to the combined Tevatron 
MI reach. In this case, there is a large Tevatron reach due to W\Z 2 — > 3£ extending to 
mi/2 ~ 230 GeV, overlapping considerably with the MDM region. Finally, we note once 
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again that LHC can cover the whole plane via multi-jet+^ T searches. In addition, the 
MDM region lies again almost entirely within the LHC slepton search region, and overlaps 
substantially with the LHC WiZ 2 — > 3£ clean trilepton region . 

Last of all, we turn to Fig. 9, which compares the neutralino relic density with collider 
search regions for large tan/3 = 10, with [i < 0. In this case, we note that the MDM favored 
region lies entirely above the region that is searchable at LEP2. In addition, for this case, 
the lightest Higgs boson has mass rrih ^ 90 GeV throughout the plane, beyond the reach of 
LEP2 at -y/i = 190 GeV. Hence, if tan (3 is large, and the MDM scenario is correct, then there 
would be little hope of seeing SUSY at LEP2. In this case, the prospect for minimal SUGRA 
at Tevatron MI is even worse, except for the narrow region extending along m ~ 100 GeV, 
which enters into the cosmologically favored MDM region. Finally, we note that once again 
the LHC slepton reach contour excloses most of the MDM region, with the main exception 
being the band of allowed MDM region extending to large mo along mi/2 ~ 160 — 170 GeV. 
The LHC W\Z 2 — > 3£ region encloses pieces of the MDM region, but leaves significant areas 
uncovered pi]. 



V. CONCLUSION 

In this paper, working within the minimal supergravity model with radiative electroweak 
symmetry breaking and universal GUT scale soft supersymmetry breaking terms, we have 
evaluated the cosmological relic density from neutralinos produced in the early universe. 
Our technique was to evaluate all lowest order neutralino annihilation Feynman diagrams 
as complex helicity amplitudes. We then performed the necessary integrations numerically, 
preserving relativistic covariance, and avoiding the usual expansion as a power series in 
velocity. While this approach might be regarded as a brute force numerical calculation, 
it does include relativistic thermal averaging, annihilation threshold effects, and careful 
integration over Breit-Wigner poles. We do not include co-annihilation processes in our 
calculations, which are however unimportant within the SUGRA framework, in which we 
work. 

Our numerical results for the neutralino relic density were presented in Figures 1-6. For 
the favored mixed dark matter scenario, for which 0.15 < Qh 2 < 0.4, we find that, unless 
annihilation occurs via s-channel Z or h exchange (in which case rrig < 300 — 400 GeV [j37|1 ), 
irtg < 1000 GeV, and < 250 GeV. The less conservative constraint from the age of the 
universe (Qh 2 < 1) yields larger bounds on sparticle masses. 

We also examined the implications of our relic density calculations for collider searches 
for the sparticles of minimal SUGRA. These results have been summarized in Figs. 7-9. 
Within the MDM range of Qh 2 , we find that LEP2 has a high probability to detect a light 
Higgs boson if tan/3 is small and /x < 0. For the opposite sign of //, m h can be larger, and 
detection at LEP2 is less certain. Prospects for detection of sleptons or charginos at LEP2 
are less bright: generally, if < 90 GeV, t-channel neutralino annihilation is too large, 
leading to rather low values of neutralino relic density. Likewise, if < 90 GeV, then 
m Zi ~ ^5 GeV, and neutralinos can annihilate via s-channel Z or h exchange, again leading 
to only a small relic abundance. 

Prospects for discovering SUGRA at Tevatron MI experiments are somewhat brighter, 
since a reasonable relic density requires roughly 100 < < 250 GeV. Such a slepton mass 
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range generally leads to enhanced leptonic decays of neutralinos, giving Tevatron experi- 
ments a good chance to find SUGRA via W\Z 2 — > 3£ searches. 

The CERN LHC pp collider can make a thorough search for supersymmetry over all the 
allowed parameter space in the multi-jet +If T channel. However, the rather light slepton 
masses required for reasonable neutralino relic densities falls within the range of LHC ex- 
perimental sensitivity, so there is a good chance to find sleptons at LHC if, for instance, 
the MDM scenario turns out to be correct. Likewise, experiments at an e + e~ linear collider 
operating at s/s ~ 500 GeV would stand a good chance of discovering sleptons, since they 
would be sensitive to slepton pair production for m„- < 230 GeV PSL 



Note added: Upon completion of this work, a preprint appeared which addressed the neu- 
tralino dark matter relic density in SUGRA models with non-universal soft-breaking terms 

US- 
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TABLES 



TABLE I. A tabulation of Feynman diagrams contributing to our neutralino annihilation cross 
section calculation. 



Process 


Particles exchanged 


s-channel 


t-channel 


u-channel 


Z1Z1 -> 


/i,ff 


^1,2,3,4 


-^1,2,3,4 


Z X Z X -> W+W 


/i, ff, Z° 


W± 2 




Z1Z1 -> Z°/i 




Z\ : 2,3,4 


^1,2,3,4 


Z1Z1 -> Z°ff 


Z°,A 


-^1, 2,3,4 


^1,2,3,4 


Z1Z1 -> 


h,H 


^1,2,3,4 


^1,2,3,4 


^iZi -> W~ff+ 


h,H,A 






^iZi -> IL+ff- 


h, ff, A 






Z1Z1 -> hh 


h : H 


^1,2,3,4 


^1,2,3,4 


Z X Z X -> ffff 


h,H 


^1,2,3,4 


^1,2,3,4 


Z1Z1 -> /iff 


h,H 


^1,2,3,4 


^1,2,3,4 


^iZi -> A A 


h,H 


-^1, 2,3,4 


^1,2,3,4 


Z1Z1 -> A 


Z\A 


Zl, 2,3,4 


^1,2,3,4 


^iZi -> HA 


h,H 


-^1, 2,3,4 


^1,2,3,4 


Z X Z X -> ff+ff- 


h, ff, Z° 








Z°, h, ff, A 
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FIGURES 



FIG. 1. Plot of contours of constant £lh 2 in the mo vs. mi/ 2 plane, where Aq = 0, tan/3 = 2, 
H < and vrit = 170 GeV. The regions labelled by TH (EX) are excluded by theoretical (experi- 
mental) considerations. 

FIG. 2. Relativistic thermally averaged cross section times velocity as a function of m^j for 
rriQ = 200 GeV, with other parameters as in Fig. 1. 

FIG. 3. Same as Fig. 1, except now \i > 0. 

FIG. 4. Relativistic thermally averaged cross section times velocity as a function of m 1( / 2 , for 
rriQ = 200 GeV, with other parameters as in Fig. 3. 

FIG. 5. Same as Fig. 1, except now tan/? = 10. 

FIG. 6. Relativistic thermally averaged cross section times velocity as a function of m^, for 
mo = 200 GeV, with other parameters as in Fig. 5. 

FIG. 7. Same as Fig. 1, except we also show the region below the contour labelled MI, which 
is accessible to Tevatron Main Injector era experiments, and the region below the LEP2 contour, 
where sparticles are accessible to experiments at LEP2 operating at \/s = 190 GeV. Below the 
LEP2-Higgs contour, the lightest SUSY Higgs h is accessible at LEP2. Below the LHC contour, 
sleptons ought to be accessible to LHC experiments. 

FIG. 8. Same as Fig. 7, except now \i > 0. 

FIG. 9. Same as Fig. 7, except now tan/? = 10. 
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